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. . ABSTRACT 

We systematically explore the evolution of the merger of two carbon-oxygen (CO) white 
i— i \ dwarfs. The dynamical evolution of a 0.9 M +0.6 M CO white dwarf merger is followed 

4^ ■ by a three-dimensional SPH simulation. The calculation uses a state-of-the art equation of 

Qh state that is coupled to an efficient nuclear reaction network that accurately approximates all 

q \ stages from helium burning up to nuclear statistical equilibrium. We use an elaborate pre- 

■ scription in which artificial viscosity is essentially absent, unless a shock is detected, and a 
^ 1 much larger number of SPH particles than earlier calculations. Based on this simulation, we 

. suggest that the central region of the merger remnant can, once it has reached quasi-static 

equilibrium, be approximated as a differentially rotating CO star, which consists of a slowly 

| rotating cold core and a rapidly rotating hot envelope surrounded by a centrifugally supported 

disc. We construct a model of the CO remnant that mimics the results of the SPH simulation 

■ using a one-dimensional hydrodynamic stellar evolution code and then follow its secular evo- 
0^ \ lution, where we include the effects of rotation on the stellar structure and the transport of 
{N| ■ angular momentum. The influence of the Keplerian disc is implicitly treated by considering 

' mass accretion from the disc onto the hot envelope. The stellar evolution models indicate that 

^-jl | the growth of the cold core is controlled by neutrino cooling at the interface between the core 

■ and the hot envelope, and that carbon ignition in the envelope can be avoided despite high 
| effective accretion rates. This result suggests that the assumption of forced accretion of cold 
. matter that was adopted in previous studies of the evolution of double CO white dwarf merger 

remnants may not be appropriate. Specifically we find that off-center carbon ignition, which 
. £h ! would eventually lead to the collapse of the remnant to a neutron star, can be avoided if the 

. following conditions are satisfied: (1) when the merger remnant reaches quasi-static equilib- 

1 rium, the local maximum temperature at the interface between the core and the envelope must 

be lower than the critical limit for carbon-ignition. (2) Angular-momentum loss from the cen- 
tral merger remnant should not occur on a time scale shorter than the local neutrino cooling 
time scale at the interface. (3) The mass-accretion rate from the centrifugally supported disc 
must be sufficiently low (M < 5 x 10 -6 ...10 -5 M yr _1 ). Our results imply that at least 
some products of double CO white dwarfs merger may be considered good candidates for 
the progenitors of Type la supernovae. In this case, the characteristic time delay between the 
initial dynamical merger and the eventual explosion would be ~ 10 5 yr. 
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1 INTRODUCTION 

The coalescence of two carbon-oxygen (CO) white dwarfs with 
a combined mass in excess of the Chandrasekhar limit has long 
be en considered a p romising p ath towards a Type la supernova (SN 
la; Iben & Tutukov 1984; Webbink 1984). Indeed, in the last few 
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years, a few massive double CO white dwarf systems have been 
found that have per iods short enough for them to merge within 
a Hubble time (e.g. iNaoiwotzki et allEx)! l2004h . This double- 
degenerate (DD) scenario can also easily explain the lack of hy- 
drogen and helium lines in most SN la spectra and the occur- 
rence of SNe la bo th in old and young star-forming systems (e.g. 
iBranch etal.ll 19951) . 

Theoretically, the final fate of double CO white dwarf merg- 
ers has been much debated. Previous studies assumed that the dy- 
namical disruption of the Roche-lobe filling secondary should lead 
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Figure 1. Schematic illustration of the configuration of the remnant of a 
double CO white dwarf merger once quasi-static equilibrium has been es- 
tablished. 



to the formation of a thick disc around the primary white dwar f 
(iTutukov & Yungelsonl 1 19791 : iMochkovitch & Livio([l98^ . Il990h . 
Therefore, accretion of CO-rich matter from the thick disc onto 
the central cold white dwarf has been studied for i nvestigating the 
evolution of such mergers by many authors jNomoto & Ibenll 19851 ; 
ISaio & Nom^ll98lTl998Ll2004l : IPiersanti et al.ll2QQ3al lbh. As ac- 
cretion rates from the thick disc should be close to the Edding- 
ton limit (M w 10 -5 M yr _1 ), most of those studies concluded 
that carbon ignition in the envelope of the accreting white dwarf 
is an inevitable consequence of such rapid accretion of CO-rich 
matter. Once carbon ignites off-center, the burning flame propa- 
gates inwards on a relatively short time scale (~ 5000 yr), and 
the CO white dwarf is trans formed into an ONeMg white dwarf 
(ISaio & Nom oto 1985. 1l998» . When the mass of the ONeMg white 
dwarf approaches the Chandrasekhar limit, electron capture onto 
Ne and Mg is expected to lead to the gravitational collaps e of 
the white dwarf to a neutron star dNom oto & Kondo 1991; see 
|PessartetaT]l2006l and lKitaura. Janka & HillebrandJ l2006l for re- 
cent studies of such collapse). 

However, the evolution of the remnants of double CO white 
dwarf mergers is not yet well understood. For instance, it has 
been debated whether the accretion rate d ecreases when the ac- 
creting white dwarf re aches critical rotation jPiersanti et al . 2003a; 
ISaio & Nomotol [20041) . More importantly, the canonical descrip- 
tion of the merger remnant as a primary white dwarf + thick 
disc system is clearly an oversimplification. In previous three- 
dime nsional smoothed particle hydrodynamics (SPH) simula- 
tions (iBenz et al.ll 19901 ; ISegretain. Chabrier & Mochkovitchlll997l ; 
iGuerrero et al.ll2004l ; see also Sect. [2]), a large fraction of the dis- 
rupted secondary and the outermost layers of the primary form an 
extended hot envelope around the cold core containing most of the 
primary mass. The rest of the secondary mass becomes a centrifu- 
gally supported disc in the outermost layers of the merger rem- 
nant. Interestingly, the merger remnant reaches a state of quasi- 
static equilibrium within a few minutes from the onset of the dy- 
namical disruption of the secondary. As the structure of the cold 
core plus the hot envelope appears to have a fairly spheroidal shape 
(see below) rather than the t oroidal shape obtained with a zero - 
temperature equation of state ( Mochkovitch & Livid 1 19891. 1 19901) . 
the merger remnant may be better described as a differentially ro- 
tating single CO star consisting of a slowly rotating cold core and 
a rapidly rotating hot extended envelope surrounded by a Keple- 
rian disc, as illustrated in Fig. [TJ than the previously adopted pri- 
mary white dwarf + thick disc system. The further evolution of the 
merger must therefore be determined by the thermal cooling of the 
hot envelope and the redistribution of the angular momentum in- 
side the central remnant, and accretion of matter onto the envelope 
from the Keplerian disc. 

With this new approach to the problem in mind, we here re- 



visit both the dynamical and the secular evolution of double CO 
white dwarf mergers. In the following section (Sect. [3, we present 
the numerical results of an SPH simulation of the dynamical evolu- 
tion of the coalescence of a 0.9 M WD and a 0.6 M CO white 
dwarf up to the stage of quasi-hydrostatic equilibrium, and we care- 
fully investigate the structure of the merger remnant. In Sect. [3] we 
construct models of the central remnant in quasi- static equilibrium 
state (primary + hot extended envelope) which mimic the SPH re- 
sult and calculate the thermal evolution of the merger remnant using 
a hydrodynamic stellar evolution code. In particular, the conditions 
for avoiding off-center carbon ignition are systematically explored. 
In Sect.lU we conclude this work by discussing uncertainties in our 
assumptions, the implications for Type la supernovae and future 
work. 



2 DYNAMICAL EVOLUTION OF THE MERGER 

Before discussing the subsequent thermal evolution after the co- 
alescence of a double CO white dwarf coalescence binary, we 
investigate the configuration of the remnant in quasi- static equi- 
librium in some detail. For this purpose, we have carried out 
a SPH simulation of the dynamical process of the coalescence 
of two CO white dwarfs of 0.9 M and 0.6 M , respec- 
tively. Our simulation uses a 3D smoothed particle hydrodynam- 
ics (SPH) code that is a n offspring of a code developed to simu- 
late neutron star mergers ([Rosswog et al.l2000l ; lRosswog & Daviesl 
120021 ; iRosswog & Lieb endorfer l2003h . It uses an artificial viscos- 
ity scheme with time-dependent parameters (M orris & Mon aghan 
Il9971) . In the absence of shocks, the viscosity parameters have a 
very low value (a = 0.05 and (3 — 0.1; most SPH implemen- 
tations use values of a — 1...1.5 and 8 — 2... 3); if a shock is 
detected, a source term dRosswog et al.l l2000) guarantees that the 
parameters rise to values that are able to resolve the shock prop- 
erly without spurious post-shock oscillations. To suppress artificial 
viscosity forces in pure shear flows, we additionally apply a switch 
originally suggested byEalsara(1995|). 

To account for the energetic feedback onto the fluid from 
nuclear transmutati ons, we use a m inimal nuclear reaction net- 
work developed bv lHix et"aD (1998). It couples a conventional a- 
network stretching from He to Si with a quasi-equilibrium-reduced 
a-network. Although a set of only seven nuclear species is used, 
this network reproduces the energy generation of all burning stages 
from H e-burning t o NSE very accurately. For details and tests we 
refer to iHix et al.l dl998h . We use the HELMHOLTZ equation of 
state (EOS), developed by the Center for Astrophysical Thermonu- 
clear Flashes at the University of Chicago. This EOS allows to 
freely specify the chemical composition of the gas and can be cou- 
pled to nuclear reaction networks. The electron/positron equation 
of state has been calculated without approximations, i.e. it makes 
no assumptions about the degree of degeneracy or relativity; the 
exact expressions are integrated numerically to machine precision. 
The nuclei in the gas are treated as a Maxwell-Boltzmann gas, 
the photons as blackbody radiation. The EOS is used in tabular 
form with densities ranging from 10 -10 ^ pY e ^ 10 11 g cm -3 
and temperatures from 10 4 to 10 11 K. A sophisticated, biquin- 
tic Hermite polynomial interpolation is used to enforce thermody- 
nami c consistency (i.e. the M axwell-relations) at interpolated val- 
ues dTimmes & S westv 20jxj). 

We use a M acCormack predictor-corrector method (e.g. 
lLomax et~al]|200ll) with individual particle time steps to evolve 
the fluid. With our standard parameters for the tree-opening cri- 
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Figure 2. Dynamical evolution of the coalescence of a 0.6 Mq + 0.9 Mq CO white dwarf binary. The panels in the left column show the density in the 
orbital plane, the panels in the right column the temperature in units of 10 6 K. Lengths are in code units (= 10 9 cm). 
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Figure 3. Dynamical evolution of the coalescence ofaO.6M0 +0.9Mq CO white dwarf binary. Continued from Fig. [2 
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Figure 4. The evolution of the local peak of temperature during the merger 
of two CO white dwarfs of 0.6 Mq and 0.9 Mq, respectively, as a function 
of time after the onset of the simulation. 



terion and the integration, this time marching implementation con- 
serves the total energy to better than 4 x 10 -3 and the total an- 
gular momentum to better than 2 x 10 -4 . Note that this could, 
in principle, be improved even further by taking into account the 
so-cal led "g rad-h"-terms dSpringel & Hernquistl 12002: Monaghan 
120021 : |Pricdl2004h and extra-terms arising from adapting gravita- 
tional smoothing terms dPrice & Monagha n 2003). 

To avoid numerical artifacts, we only use equal mass SPH 
particles. For the initial conditions, we therefore stretch a uniform 
particle distribution according to a function that has been derived 
from solving the ID stellar structure equations. This techni que is 
described in detail in lRosswog. Ramirez-Ruiz & Hixl J2007). This 
particle s etup is then further relaxed with an additional damping 
force (e.g. lRosswog. Speith & Wvnnll2004l) so that the particles can 
settle into their true equilibrium configuration. The calculations are 
performed with 2 x 10 5 SPH particles, a much larger particle num- 
ber than could be afforded by previous calculations, and run up to 
a much longer evolutionary time (5 minutes) than previous calcu- 
lations (see Table Q]). 

Figs. [2] and [3] show the dynamical evolution of the merging 
process of the double white dwarf system considered in this study. 
The panel in the left columns show the densities and the panel in 
the right columns the temperatures (in units of 10 6 K) in the orbital 
plane. The secondary is completely disrupted within 1.7 minutes, 
and mass accretion onto the primary induces local heating near the 
surface of the primary. Fig. |4] shows the evolution of the maximum 
temperature as a function of time. The peak in temperature reaches 
1.7xl0 9 Katt ~ 1.0 min, where t = 0.0 marks the moment when 
the simulation starts. Carbon ignites when T > 10 9 K, but nuclear 
burning is quenched soon due to the local expansio n of the hottest 
region, as is also observed in the simulations of iGuerrero et al.l 
(2004). The total amount of energy released due to nuclear burn- 
ing i s about 10 45 erg. 

ISegretain. Chabrier & Mochkovitchl dl997h considered the 
same initial white dwarf masses as in the present study. But 
they adopted the origina l artificial viscosity prescription of 
iMonaghan & Varnasl dl988l) . which is known to introduce spuri- 
ous forces in shear flows, and they did not include nuclear burning 
(Table [TJ- By the end of their calculation (t — 1.56 min), T max 
reached 8 x 10 8 K, while in our simulation, it decreases to 8 x 10 8 K 
only when £ ^ 1.7 min. Interestingly, T max decreases further af- 
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Figure 5. Top: Density contour of the merger remnant in the x — z plane at 
t = 5.3 min. Here one code unit corresponds to 10 9 cm. Middle: Thermo- 
dynamic structure of the merger remnant at t = 5.31 min: shown are the 
temperature and the density as a function of distance from the centre, along 
the positive x- and z-axis, as indicated. Bottom: Angular velocity in units 
of the local Keplerian value at t = 5.31 min, along the positive/negative 
x- and y-axis of the merger remnant. 



terwards in our calculation, as shown in Fig.|4] and reaches a steady 
value at T ma x ~ 5.6 x 10 8 K whe n t > 2.5min. In the o ther cal- 
culations bv lBenz etail Jl99Qh and lGuerrero et all (l2004h . the dy- 
namical evolution of the merger was not followed for more than 2 
minutes either, and we cannot directly compare our results to theirs. 
However, we suspect that T max would also decrease further in the 
systems they considered if they had continued their calculations for 
a longer evolutionary time. It should also be noted that energy dis- 
sipation by artificial viscosity might lead to overheating, and that 
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thermal diffusion - which may play an important role in the outer- 
most layers - is not considered in the present study. It is thus likely 
that T m ax in the quasi- static equilibrium state may be even lower 
in reality than in our simulation. 

Fig.[5]shows the structure of the merger remnant at quasi-static 
equilibrium. The central region with R < 10 9 cm (M r < 1.1 M©) 
has a fairly spheroidal shape, and a centrifugally supported disc 
appears at R > 10 9 cm where the angular velocity is close to 
the Keplerian value. The fraction of the secondary mass contained 
in the Keplerian disc is larger in our simulati on (about 67%) than 
in ISegretain. Chabrier & Mochkovitchl dl997h (about 41 %). The 
innermost core (R < 3 x 10 8 cm; M r < 0.6 M©) is essen- 
tially isothermal, and the temperature has its peak value (T ~ 
5.6 x 10 8 K) at R ~ 5 x 10 8 cm and M r ~ 0.85 M©. The 
disc material extends over 4 x 10 9 cm along the z-axis as the tem- 
perature is still high; if thermal diffusion were included, the disc 
would become much thinner on a short time scale of a few hours. 
Therefore, our simulation confirms the remnant structure at quasi- 
static equilibrium that is illustrated in Fig.Q] In the next section, we 
investigate the secular evolution of the merger from such a quasi- 
static equilibrium state. 

3 SECULAR EVOLUTION OF THE MERGER REMNANT 
3.1 Physical assumptions and methods 

Our SPH simulation shows that the remnant of the merger of two 
CO white dwarfs (0.9 M© + 0.6 M©) in the state of quasi-static 
equilibrium has the following features (see Fig. [6]): 

(i) The core is cold and nearly isothermal. 

(ii) The local peak of temperature (T p ) is located at a mass co- 
ordinate slightly less than the primary mass. 

(iii) A steep gradient in temperature appears at the interface be- 
tween the core and the local peak of temperature. 

(iv) The interface is rather widely extended into the primary 
(AMinterface ~ 33 % of the primary mass), and the mass of the 
quasi-isothermal cold core (M CO re) is about 77 % of the primary 
mass. 

(v) The mass of the outer envelope above the local peak of tem- 
perature contains about 33 % of the mass of the secondary, and the 
rest of the secondary forms a Keplerian disc. 

Let us define T p as the local peak of temperature at quasi- static 
equilibrium, Mcm as the mass of the central remnant (cold core + 
hot envelope), and M p as the location of T p in the mass coordi- 
nate (i.e., M p = Mcore + AMinterface; see Fig. O. To construct 
models of the central remnant, we use a one-dimensional hydrody- 
namic stellar evolution code which incorporates the effects of ro- 
tation on the stellar structure, transport of angular momentum due 
to the shear instability, Eddington- Sweet circulation, and the Gol- 
dreich, Schubert and Fricke instability, and dissipation of rotational 
energy due to shear motions. The effects of magnetic fields are ne- 
glected (see Se c. [4]). More details about the code are described in 
(lYoon & Langel2004l : hereafter YL04) and references therein. 

In order to mimic the temperature profile of the central rem- 
nant as obtained from the SPH simulation, we artificially deposit 
energy in the envelope, using the following prescription for a white 
dwarf with M — Mcr : 

e (M r ) = A(T'(M r ) - T(M T )) [erg g" 1 s" 1 ], (1) 
where 



- 1 1 ' 






j 


/T\ 


j 


= 


/ M o = \ 


= 




/ 0.84M G \ 






M core =0.6M Q / 






^ \J < AM env > 




= 1 1 1 


N M interf"5 




0.0 


0.2 0.4 0.6 0.8 1.0 






M r /M Q 






5 10 15 

R [10 8 cm] 

0.6 1 1 1 I 1 1 1 1 I 1 1-0 




o.i E_. i i i i i i i i i i i i i i i Jo.o 

5 10 15 

R [10 8 cm] 



Figure 6. Initial model of the central remnant for sequences Sal - Sail. 
The top and middle panels show temperature as a function of the mass co- 
ordinate and radius, respectively. The solid curve in the bottom panel gives 
the angular- velocity profile as a function of radius. The dashed curve de- 
notes the angular velocity in units of the local Keplerian value. 



T'(M r ) 

3-10 7 K+(7-10 7 K-3-10 7 K)( I ^ r ) 2 , 

ifM r < M C ore, 

< T p -(T p -7.10 7 K)( ] ^ r ) 2 , (2) 

if M core < M r < M p , 

r — (r — T \ ( iog[p(M r )/p s ] \ 2 

° V° J-P) \lo g [p(M p )/ Ps ]J ' 

if M r > M p . 

In this way, the temperature profile in the central remnant model 
follows T f (M T ). Here, A and C are constants. We use A = 
10 5 ergg -1 s" 1 K" 1 and C = 2 x 10 8 K in most cases. 
A rotational profile is imposed as 
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Table 1. Comparison of SPH simulations of double CO white dwarf mergers. The columns list: Mi and M 2 : the masses of t he primary and the secondary, 
respectively; NoP: the total number of particles used; ^ sp h: the type of artificial viscosity employed, 'std.' refers to Monaghan & Varnas (1988); Network: 
type of nuclear network employed; t s i m : evolutionary time that has elapsed by the end of the calculation; T max : maximum temperature obtained during the 
simulation; and T p : the local peak of temperature at the end of the calculation 



Ref.* 


Mi 


M 2 


NoP 


^sph 






Network 


^sim 


Tmax 


T P 




1 


1.2 M 


0.9 M 


~ 7 x 10 3 


std. 








None 


51 sec. 


7 


~ 10 9 K 


2 


0.8 M 


0.6 M 


~ 4 x 10 4 


std. 


+ Balsara-switch 






alpha network 


50 sec. 


1.4 x 10 9 K 


? 




2 


l.OM 


0.6 M 


~ 4 x 10 4 


std. 


+ Balsara-switch 






alpha network 


65 sec. 


1.6 x 10 9 K 


? 




2 


1.0 M© 


0.8 M 


~ 4 x 10 4 


std. 


+ Balsara-switch 






alpha network 


65 sec. 


2.0 x 10 9 K 


? 




3 


0.9 M 


0.6 M 


- 6 x 10 4 


std 








None 


1.56 min. 


? 


~ 7 x 10 8 K 


4 


0.9 M 


0.6 M 


2 x 10 5 


see 


Rosswos et al. (2000) 


QSE-alpha network 5.3 min. 


1.7 x 10 9 K 


5.6 x 10 8 K 
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Table 2. Merger remnant model sequences. Each column lists the following: 


No 


: sequence label, Mcr: mass of the central remnant, M CO re: 


mass of the 


quasi-isothermal core 


M p : location of the local peak of temperature in the mass coordinate, T p 


the local peak of temperature, p p : density at M r = M p , 


tj: adopted time scale for angular momentum loss according to Eq. (4), M acc 


: adopted mass accretion rate from the Keplerian disc, C-ig: off-center ignition 


of carbon, Mwd in- 


total mass of the central remnant when off-center carbon 


ignition occurs, 


M r i g : location of off-center carbon ignition 


in the mass 


coordinate 


























No. 


M CR 


-Mcore 


M p 




T P 


Pp 




rj 


Macc 


C-ig 


A^WD,ig 


M r;ig 




[M ] 


[M ] 


[M ] 




[10 8 K] [10 6 g/cm c 


5 ] [yr] 


10- 6 M /yr 




[M ] 


[M ] 


Sal 


1.11 


0.6 


0.84 




5.6 


0.8 




00 


0.0 


No 


- 


- 


Sa2 


1.11 


0.6 


0.84 




5.6 


0.8 




10 2 


0.0 


Yes 


1.11 


0.80 


Sa3 


1.11 


0.6 


0.84 




5.6 


0.8 




10 3 


0.0 


Yes 


1.11 


0.80 


Sa4 


1.11 


0.6 


0.84 




5.6 


0.8 




10 4 


0.0 


Yes 


1.11 


0.85 


Sa5 


1.11 


0.6 


0.84 




5.6 


0.8 




10 5 


0.0 


No 


- 


- 


Sa6 


1.11 


0.6 


0.84 




5.6 


0.8 




10 5 


10.0 


Yes 


1.34 


1.09 


Sa7 


1.11 


0.6 


0.84 




5 6 


8 




10 5 


5.0 


Yes 


1.34 


1.20 


Sa8 


1.11 


0.6 


0.84 




5.6 


0.8 




10 5 


2.0 


No 


- 


- 


Sa9 


1.11 


0.6 


0.84 




5.6 


0.8 




10 5 


1.0 


No 


- 


- 


SalO 


1.11 


0.6 


0.84 




5.6 


0.8 




5- 10 5 


5.0 


No 




- 


Sail 


1.11 


0.6 


0.84 




5.6 


0.8 




5 • 10 5 


1.0 


No 






Aal 


1.25 


0.6 


0.93 




5.0 


2.3 




00 


0.0 


No 






Aa2 


1.25 


0.6 


0.93 




5.0 


2.3 




10 2 


0.0 


Yes 


1.250 


0.90 


Aa3 


1.25 


0.6 


0.93 




5.0 


2.3 




10 3 


0.0 


Yes 


1.250 


0.92 


Aa4 


1.25 


0.6 


0.93 




5.0 


2.3 




10 4 


0.0 


Yes 


1.250 


1.12 


Aa5 


1.25 


0.6 


0.92 




5.0 


2.3 




10 5 


0.0 


No 






Aa6 


1.25 


0.6 


0.92 




5.0 


2.3 




10 6 


0.0 


No 






Aa7 


1.25 


0.6 


0.92 




5.0 


2.3 




10 5 


10.0 


Yes 


1.360 


1.20 


Aa8 


1.25 


0.6 


0.92 




5.0 


2.3 




10 5 


5.0 


No 






Aa9 


1.25 


0.6 


0.92 




5.0 


2.3 




10 5 


1.0 


No 






AalO 


1.25 


0.6 


0.92 




5.0 


2.3 




10 6 


10.0 


Yes 


1.382 


1.22 


Abl 


1.25 


0.7 


0.92 




5.0 


3.1 




00 


0.0 


No 






Ab2 


1.25 


0.7 


0.92 




5.0 


3.1 




10 3 


0.0 


Yes 


1.250 


0.97 


Ab3 


1.25 


0.7 


0.92 




5.0 


3.1 




10 4 


0.0 


No 






Ab4 


1.25 


0.7 


0.92 




5.0 


3.1 




10 5 


0.0 


No 






Ab5 


1.25 


0.7 


0.92 




5.0 


3.1 




10 5 


10.0 


Yes 


1.344 


1.21 


Ab6 


1.25 


0.7 


0.92 




5.0 


3.1 




10 5 


5.0 


No 






Acl 


1.25 


0.5 


0.88 




5.9 


1.6 




00 


0.0 


Yes 


1.250 


0.84 


Adl 


1.25 


0.6 


0.92 




6.0 


1.8 




00 


0.0 


Yes 


1.250 


0.90 


Ad2 


1.25 


0.6 


0.92 




6.0 


1.8 




10 6 


5.0 


Yes 


1.252 


0.90 


Ael 


1.25 


0.6 


0.90 




6.8 


1.5 




00 


0.0 


Yes 


1.250 


0.87 


Bal 


1.363 


0.82 


0.95 




5.0 


12.2 




00 


0.0 


No 






Ba2 


1.363 


0.82 


0.95 




5.0 


12.2 




10 2 


0.0 


Yes 


1.363 


0.95 


Ba3 


1.363 


0.82 


0.95 




5.0 


12.2 




10 3 


0.0 


Yes 


1.363 


1.12 


Ba4 


1.363 


0.82 


0.95 




5.0 


12.2 




10 4 


0.0 


No 






Ba5 


1.363 


0.82 


0.95 




5.0 


12.2 




10 5 


0.0 


No 






Ba6 


1.363 


0.82 


0.95 




5.0 


12.2 




10 5 


10.0 


Yes 


1.398 


1.34 


Ba7 


1.363 


0.82 


0.95 




5.0 


12.2 




10 5 


5.0 


Yes 


1.483 


1.43 


Ba8 


1.363 


0.82 


0.95 




5.0 


12.2 




10 5 


1.0 


No 






Tal 


1.25 


0.60 


0.86 




5.0 


28.8 




00 


0.0 


No 







8 S.-C. Yoon, Ph. Podsiadlowski & S. Rosswog 




234567 23456 
log p [cm g" 3 ] log p [cm g" 3 ] 



Figure 7. (a) Evolution of a non-rotating white dwarf accreting with a constant accretion rate of M = 10 _ 5 M0 yr _1 with an initial mass of 0.9 Mq 
(Seq. N0.9) in the density - temperature plane, (b) The local effective accretion rate (M e ff ,r : = 47rr 2 pv) as a function of density in Seq. N0.9, at different 
evolutionary epochs as indicated by the labels, (c) - (f) The rates of energy loss/production due to neutrino (e u ) cooling, compressional heating (e com p), 
nuclear energy generation (e nuc ) and thermal diffusion (e t h) at different evolutionary epochs. Note that here e u , e CO mp and e nU c represent the values which 
are used in the evolutionary calculations, while e t h is an order-of-magnitude estimate according to Eq. (7). 



Q(M r ) 

c^o + (1 



)( 



M r -M cor , 
M CK -M co 



if M r < M core , 
GM r 

if M r ^ M C ore, 



(3) 



where Q = 0.2^GM C r/R 3 and u = Ho/ \/ GM cove /rl ore . 
As shown in Fig. [6] this simple assumption gives a rotational ve- 
locity profile that is morphologically similar to that found in the 
SPH simulation: a steep gradient at the interface between the core 
and the envelope, and a local peak in the envelope. Within our 1-D 
approximation of the effects of rotation, the exact shape of the ro- 
tational velocity profile does not affect the main conclusions of the 
present work for the following reasons. Firstly, the velocity gradient 
at the interface is adjusted to the threshold value for the dynamical 
shear instability on a very short time scale (see below, and discus- 



sions in YL04). Secondly, our 1-D approximation underestimates 
the effect of the centrifugal force on the stellar structure in layers 
which rotate more rapidly than about 60 % critical (YL04), and un- 
certainties due to this limit are much greater than due to the shape 
of Q(r) in the outer layers of the envelope. Possible uncertainties 
due to this limitation are critically discussed in Sect.|4] 



The central remnant may lose angular momentum by out- 
ward angular momentum transport into the Keplerian disc 
JPopham & Naravanll 199 ll : IPaczvnskill 199 lh and/or by the gravita- 
tional wave radiation, e.g., due to th e r-mode instability (I Anders son! 
1998; Friedman & Morsink 1998). Our code cannot properly de- 
scribe any of these effects, and here we consider them simply by 
ass uming a const ant time scale for the ang ular momentum loss (tj; 
see Knaapl l2004l : cf. IPiersanti et al . 2003a), such that the specific 
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Table 3. Accreting white dwarf model sequences with a constant accretion 
rate (M = 10 -5 M /yr). The columns list: No: sequence label, Mi n i t : 
initial mass, logLi n i t /L : initial luminosity, MwD,ig : the total mass of 
the white dwarf when carbon ignites off-center, and M r> i g : location of car- 
bon ignition in the mass coordinate. Sequences with 'N' are for non-rotating 
models, and 'R' for rotating models. 
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N0.9 
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Figure 8. Local effective mass accretion rate (M e ff,r = 4irr 2 pv) as a 
function of density in the models of sequence R0.9, at different evolutionary 
epochs. 



angular momentum of each mass shell decreases over a time step 
At by an amount 

Aji — ji [1 — exp(— At/rj)] . (4) 

Mass accretion from the Keplerian disc is also considered in 
some model sequences, with different values for the accretion rate 
(M acc ). The angular-momentum accretion is treated in the same 
way as in YL04: the accreted matter is assumed to carry angular 
momentum at a value close to the Keplerian value if the surface 
velocity of the central remnant is below critical, while no angular- 
momentum accretion is allowed otherwise. 

Model sequences with different sets of Mcr, M CO re, M p , rj, 
M, and T p are calculated, as summarized in Table The initial 
model in Seqs S is intended to reproduce the result of our SPH sim- 
ulation, where M C r = 1.10 M and M p w 0.84 M are adopted. 
We also assume M C r = 1.25 M and M p & 0.9 M in Seq. A, 
and Mcr = 1.364 M and M p = 0.95 M in Seq. B, to simulate 
mergers of 0.9 - 1.0 M + 0.7 - 1.0 M white dwarf binaries. 
At a given Mcr, different sets of M core , M p , and T p are marked 
in the sequence label by minor characters (a, b, c, d, e), while dif- 
ferent sets of rj and M acc are indicated by Arabian numbers. For 
instance, sequences Sal-Sail have the same initial merger model, 
but different values for rj and M acc - Rotation is neglected in a test 
sequence Tal (i.e., the models are non-rotating). The temperature 
and angular- velocity profiles in the initial central remnant model of 
Seqs Sal - Sal 1 are shown in Fig.[6] The temperature (a few to sev- 
eral 10 8 K) and the size (~ 10 9 cm) of the envelope appear to be 
comparable to those obtained from the SPH simulation (see Fig. [5]). 

For comparison, we also ran model sequences for classi- 
cal cold-matter accretion with a constant accretion rate of M = 
10 -5 M /yr, for both non-rotating and rotating cases, as summa- 
rized in Table [3] 

3.2 Results 

3. 2. 1 Classical models of cold-matter accretion 

Before discussing the central remnant models, let us first investi- 
gate the evolution of classical cold-matter accreting white dwarf 
models in detail. In these models, the accreted matter is assumed to 
have the same entropy as the surface valu e of the accreting whit e 
dwarf. As shown in previous studies (e.g. iNomoto & Iberj|l985l) . 
the thermal evolution of rapidly accreting white dwarfs is deter- 



mined by the interplay of compressional heating and thermal dif- 
fusion. Fig. [7] shows an example of the evolution of such accret- 
ing white dwarf models for an initial WD mass of 0.9 M and 
a constant accretion rate of M acc = 10 -5 M yr -1 (Seq. N0.9; 
Tabled. 

Fig. [7^ shows that the temperature increases continuously in 
the envelope (p < 10 6 g cm -3 ), and finally carbon burning be- 
comes significant at p c± 5.6 x 10 5 g cm -3 and T~6x 10 8 K 
when t ^ 1.3 x 10 4 yr. In Figs.jTJ; -f, the rates of compressional 
heating (e C om P ), neutrino cooling (e u ), nuclear energy generation 
(e nU c) and thermal diffusion (eth) are shown. In our stellar evolu- 
tion code, the compressional heating rate is calculated according 
to 

e " - (^) (5) 

(Kippenhahn & Weigert 1990). Neutrino cooling rates are obtained 
following lltoh et al.M l996). While fcomp) tvi and e n uc in the figures 
correspond to the values that are used for the evolutionary calcula- 
tions, the thermal diffusion rate (e t h) - which is only calculated 
implicitly in the code - can only be estimated to within an order- 
of-magnitude from 

e t h ~ TCp/nh • (6) 

Here Cp denotes the specific heat at constant pressure, and r t h the 
local thermal diffusion time scale defined as 

r th = Hi IK , (7) 

where Hp is the pressure scale height, and K 
[(4acT 3 )/(3C P Kp 2 )] is the thermal diffusivity. It is clear 
from Fig. [7] that the local peak of temperature is located where the 
compressional heating rate begins to dominate over the thermal 
diffusion rate (p w 10 5 g cm -3 ), as expected. The neutrino 
cooling rate also increases as the temperature in the envelope 
becomes higher, but nuclear energy generation becomes significant 
before neutrino cooling dominates the thermal evolution, inducing 
a carbon-burning flash around p ~ 5.6 • 10 5 g cm -3 . 

As Table [3] show s, and consistent with the findings of 
INomoto & Ibenl dl985h . such off-center carbon flashes occur re- 
gardless of the initial mass of the white dwarf, if M aC c ~ 
10 -5 M yr -1 . The results with models including rotation show 
that carbon ignition may be delayed if the eff ect of rotation is in- 
clude d (Table(3j see also lPiersanti et al.ll2003al and lSaio & Nomotol 
2004). The reason is that the local effective mass accretion rate 
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Figure 9. (a) Evolution of the central remnant in Seq. Sal in the log p — T plane. The dotted curve gives the critical temperature where the nuclear energy 
generation rate due to carbon burning equals the energy loss rate due to neutrino cooling, (b) The local effective accretion rate (M e ff,r = 47rr 2 pv T ) as a 
function of density in the merger remnant model of Seq. Sal, at different evolutionary epochs as indicated by the labels, (c) - (f) The rates of energy loss/gain 
due to neutrino (e v ) cooling, compressional heating (e com p), nuclear energy generation (e comp ) and thermal diffusion (e t h) as a function of density in the 
central remnant models of Seq. Sal at different evolutionary epochs. Note that here e u , e CO mp and e nU c represent the values which are used in the evolutionary 
calculations, while e t h is an order-of-magnitude estimate according to Eq. (7). 



(M e ff ,r = 4:7rr 2 pv) inside the white dwarf at a given mass is lower 
because of the centrifugal force. For instance, in Seq. N0.9, we 
have M e ff,r ~ 10" 5 M yr -1 at around p = 5 x 10 5 g cm -3 
when t ~ 10 4 yr (Fig. [TJd), but M e ff , r is lowered by a factor of 
two in the corresponding rotating model at a similar epoch (i.e., 
M e ff,r ~ 5 x 1CT 6 M yr -1 ), as revealed in Fig. [8] However, 
carbon ignition occurs well before the white dwarf reaches the 
Chandrasekhar limit, in all model sequences considered. Thus, ro- 
tation by itself cannot change the conclusion of the previous work 
that the coalescence of double CO white dwarfs should lead to 
accretion-induced collapse rather than a thermonuclear explosion, 
un less the accretio n rate is sign ificantly lowered, as w as also shown 
bv lPiersanti etail <2003ah and lSaio & Nomotd b004h . 



3.2.2 Sequences without angular-momentum loss and mass 
accretion 

Having understood the physics of the thermal evolution of CO 
white dwarfs which accrete cold matter with a rate close to the Ed- 
dington limit, we now investigate the evolution of the central rem- 
nant model consisting of a cold core and a hot envelope as described 
in Sect. 13.11 First, we examine the results of the model sequences 
where both angular-momentum loss and mass accretion from the 
Keplerian disc are neglected (i.e., r j = oo and M — 0; Seqs Sal, 
Aal, Abl, Acl, Adl, Ael, Bal, & Tal). 

Fig. [9^ illustrates the evolution of the central remnant for 
Mcr = 1.10 M in Seq. Sal in the density - temperature 
plane. Note that the local peak of temperature at t — 0.0 (T p = 
5.6 x 10 8 K) is significantly below the critical temperature for car- 
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bon ignition (Tc-i g ; dotted curve in Fig[£k). It is shown in Fig. [9]} 
that the local effective accretion rate (M e ff,r) remains relatively 
high (5 x 10" 6 - 10" 5 M yr -1 ) around p = 10 6 g cm -3 , where 
the local peak of temperature is located, for about 5000 yrs. Despite 
such high effective accretion rates, the temperature peak continu- 
ously decreases, although the inner core becomes somewhat hotter 
due to compression, and the central remnant finally becomes a cold 
white dwarf. A few remarkable differences compared to the stan- 
dard accreting white dwarf models are found in this regard. Firstly, 
since the envelope is very hot, neutrino cooling - in particular by 
photoneutrinos - is significant from the beginning, and even domi- 
nant over the thermal diffusion at the interface between the core and 
the envelope as shown in Fig. |9j:. In cold-matter accreting white 
dwarfs, neutrino cooling becomes important only after a significant 
amount of mass has been accreted (Fig. [7]). Secondly, the compres- 
sional heating rate is slightly lower than the neutrino cooling rate 
around the local peak of temperature. As the contraction of the cen- 
tral remnant is mainly determined by the thermal evolution of the 
envelope, the local accretion rate is in fact controlled by the cool- 
ing process. This explains why we have e com p ~ e u around the 
local peak of temperature for the initial ~ 10 4 yrs, and why the 
local peak of temperature continuously decreases despite the rela- 
tively high effective accretion rate. This conclusion is the same for 
all other sequences with a Tp that is significantly lower than Tc-i g 
(Seq. Aal, Abl, & Bal) including the non-rotating case (Seq. Tal; 
Table©. 

We find that, in Seq. Sal, the differentially rotating layers at 
the interface between the core and the envelope in the initial central 
remnant model are stable against the dynamical shear instability 
(DSI). They are, however, unstable to the DSI in other sequences, 
where the interface is more degenerate (see YL04 for discussions 
on the DSI). Consequently, in Seq. Aal for example, the rate of 
rotational energy dissipation (e ro t) appears to be very high initially 
(Fig. [TO]). The differentially rotating layers are rapidly smeared out 
by the dynamical shear instability (see the discussion in Sect. 2 in 
YL04), and e ro t falls below the thermal diffusion and/or neutrino 
cooling rate only within 20 yrs. Hence we conclude that the rota- 
tional energy dissipation does not play an important role for the 
long-term evolution of the central remnant. 

Fig. [TT] shows how the evolution of the central remnant 
changes if the local peak of temperature in the initial model (T p ) is 
close to or above the critical limit for carbon burning (Tc-i g ), with 
Seqs Acl and Ael as examples. In contrast to Seq. Sal or Aal, car- 
bon burning dominates the evolution very soon in both sequences, 
and the temperature increases rapidly. Although the further evolu- 
tion has not been followed in the present study, it is most likely that 
the carbon-burning flame propagates inward such that the central 
remnant is converted into an ONeMg white dwarf within several 
thousand years as shown bv lSaio & Nomotol Jl998h . 

As summarized in Table© all other sequences follow the same 
evolutionary pattern: off-center carbon ignition is avoided in Seqs 
Abl, Bal, and Tal where T p is significantly below Tc-i g , while 
carbon ignites off-center in the other sequences where T p > Xc-i g . 
It is thus remarkable that the thermal evolution of the central rem- 
nant is sensitively determined by the local peak of temperature in 
the quasi- static equilibrium state. 

In conclusion, in the absence of angular-momentum loss and 
mass accretion from the Keplerian disc, the thermal evolution of 
the central remnant is roughly controlled by neutrino cooling at the 
interface between the core and the envelope, and off-center carbon 
burning may be avoided as long as T p < Xc-i g , while it seems 
inevitable if T p > T C -i g . 
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Figure 10. Upper panel: The angular velocity relative to the local maximum 
as a function of the mass coordinate in the central remnant models of Seq. 
Aal at different evolutionary epochs. Lower panel: the rate of rotational 
energy dissipation (e ro t ; see YL04) as a function of the mass coordinate in 
the corresponding models shown in the upper panel. 



3.3 Effect of angular momentum loss 

In Seqs Sa2 - Sa5, the central remnant has the same initial condi- 
tions as in Seq. Sal, angular momentum loss from the white dwarf 
with different time scales rj is considered according to Eq. (4). 
Note that off-center carbon ignition occurs in Seqs Sa2, Sa3 & 
Sa4, where tj < 10 4 yr, while it is avoided in Seq. Sa5 where 
rj = 10 5 yr. These results indicate that off-center carbon igni- 
tion should be induced if the angular-momentum loss occurs too 
rapidly for neutrino cooling or thermal diffusion to control the ef- 
fective mass accretion. For instance, Fig.[T2lshows that in Seq. Sa4, 
where tj = 10 4 yr, the effective mass accretion rate reaches a few 
10 -5 M yr -1 at the interface between the core and the envelope 
(p « 10 6 g cm -3 ), and the compressional heating rate exceeds the 
neutrino cooling rate. 

It is shown that the critical angular-momentum-loss time 
scale, tj, for off-center carbon ignition (rj, C rit) is smaller for Seqs 
Ab and Ba than for Seqs Sa and Aa: rj, C rit ~ 10 3 for Seqs Ab 
and Ba, and rj, C rit ~ 10 4 for Seqs Sa and Aa. This is due to the 
different local thermodynamic properties at the interface between 
the core and the envelope in different central remnant models. As 
shown in Fig.[l3] higher density and/or temperature at the interface 
result in a shorter neutrino cooling time, making it possible to avoid 
local heating for a smaller tj. In other words, Tj, cr it roughly cor- 
responds to the time scale for neutrino cooling at the local peak of 
temperature (t U)P ). 

From this experiment, we conclude that, in the absence of 
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Figure 11. Evolution of the central remnant in Seqs Acl (upper panel) and 
Ael (lower panel), in the log p — T plane. The dotted curve gives the critical 
temperature where nuclear energy generation rate due to carbon burning 
equals to energy loss rate due to neutrino cooling. 



mass accretion from the Keplerian disc, carbon ignition may be 
avoided in the central remnant, if T ma x,init < Tc-i g , and if 

Tj > T u ,p. 

3.4 Mass accretion from the Keplerian disc 

In reality, mass accretion from the Keplerian disc onto the central 
remnant is expected. The accretion rate is determined by the vis- 
cosity of the disc, which is not well known. However, we expect 
the accretion rate from a Keplerian disc may be significantly lower 
than from a pressure-supported thick disc that was assumed in pre- 
vious studies. Our results, as summarized in Tabled indicate that 
even with mass accretion, the central remnant with T p < Tc-i g can 
avoid off-center carbon ignition if the accretion rate is sufficiently 
low (i.e., M < 5 x 10~ 6 ...10~ 5 M yr" 1 ), and if rj > t u , p (see 
TableH. 

The thermal history of the central remnant in those sequences 
where carbon ignites off-center is similar to that of the white dwarf 
in classical accretion model sequences. However, as the central 
remnant has a rapidly rotating hot envelope, carbon ignition is sig- 
nificantly delayed compared to the case of classical accretion. In 
Seq. N1.2, where M in it = 1.2 M and M acc = 10~ 5 M yr" 1 , 
carbon ignites only when about 0.025 M is accreted, while in 
Seq. Aa7 more than 0.15 M have to be accreted to induce carbon 
ignition at the same accretion rate, despite its higher initial mass. 
On the other hand, the comparison of Seq. Aa7 with Seq. AalO 
indicates that off-center carbon ignition is delayed if the central 
remnant keeps more angular momentum. The critical accretion rate 



for inducing off-center carbon ignition is thus difficult to precisely 
determine, as our 1-D models significantly underestimate the effect 
of the centrifugal force, especially in the envelope where carbon 
ignites. In addition, the physics of angul ar momentum loss/g ain is 
not well understood yet, as discussed in lYoon & Lan ger (2005). 

Note that MwD,i g in Seqs AalO, Ba6 and Ba7 is already very 
close, or even above the Chandrasekhar limit. However, the cen- 
tral density in those models is still smaller by an order of mag- 
nitude than the critical limit for carbon ignition due to the effect 
of rotation. As the carbon-bur ning flame will propa gate inwards 
within several thousand years dSaio & Nomotdfl998r) . only about 
~ 0.05 M may be further accreted by the time the burning flame 
reaches the center, and the central density may not become high 
enough to induce a thermonuclear explosion before the whole cen- 
tral remnant is converted into an ONeMg white dwarf. (Super-) 
Chandrasekhar mass ONeMg white dwarfs produced in this way 
will eventually coll apse to a neutron star (see lYoon & La nger 2005; 
iDessart et ai]l2006h . 

On the other hand, the white dwarf continuously grows 
to/above the Chandrasekhar limit 1.4 M©) without suffering 
carbon ignition (neither at the center nor off-center) in Seqs Sa8, 
Sa9, SalO, Sail, Aa8, Aa9, Ab6, and Ba8. The outcome in these 
cases is thus the formation of a (super-) Chandrasekhar mass CO 
white dwarf, which will eventually explode as a Type la super- 
nova. The mass of the exploding white dwarf should d epend on 
the amount of angular momentum (lYoon & La nger 200S) and can- 
not exceed the mass budget of merging white dwarfs. Fig.[T4lshows 
the evolutionary paths of the central remnant for Seqs Sa8, Sa9 & 
Sal 1 as examples in the mass - angular momentum plane. Note that 
the central remnant initially has a large amount of angular momen- 
tum (J — 1.11 x 10 50 erg s), such that without loss/gain of angular 
momentum, it should accrete matter until it reaches M c± 1.68 M 
where it explodes in a SN la explosion. In Seqs Sa8 and Sa9, the ac- 
cretion time scale (r aC c) is longer than the angular momentum loss 
time scale, and the total angular momentum of the white dwarf con- 
tinuously decreases while the total mass increases. Consequently, 
carbon ignites at the center when the white dwarf grows to 1 .50 M 
and 1.42 M for Seqs Sa8 and Sa9, respectively. In Seq. Sail, on 
the other hand, both mass and angular momentum of the central 
remnant continuously increase, given that r aC c < tj, and a SN la 
explosion is expected only when M c± 1.70 M . Note that this 
is even larger than the mass budget of the binary system consid- 
ered for this sequence (i.e, 0.9 M + 0.6 M©). In nature, the white 
dwarf must stop growing in mass when M = 1.5 M , and a SN la 
explosion will be induced only when a sufficient amount of angular 
momentum has been removed, e.g. via gravitational wave radiation, 
as illustrated by the path Sal 1-B in Fig.[T4l 



4 CONCLUSION AND DISCUSSION 

We have explored the dynamical and secular evolution of the 
merger of double CO white dwarf binaries whose total mass ex- 
ceeds the Chandrasekhar limit. Based on our new SPH simula- 
tion of the coalescence of two CO white dwarfs of 0.9 M© and 
0.6 M©, we suggest that the immediate post-merger remnant is 
best described as a differentially rotating CO star consisting of a 
slowly rotating cold core and a rapidly rotating hot envelope that is 
surrounded by a Keplerian disc rather than as "cold white dwarf + 
thick disc" system, as in previous investigations. The evolution of 
such a CO star is determined by the thermal evolution of the en- 
velope, and the growth of the core is controlled by the cooling due 
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Figure 12. Same as in Fig. [9] but for Seq. Sa4. 
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Figure 13. Contour lines of the neutrino cooling time scale (r u = 
TCp/eu) in the logp — T plane. The level at each line gives logTv in 
units of years. The dashed line denotes the critical temperature for carbon 
ignition. The local peak of temperature and the corresponding density in the 
initial model of Seqs Sa, Aa, Ab, Ac, Ae, and Ba are marked by the filled 
symbols as indicated by the labels. 



to neutrino emission and thermal diffusion, which is fundmentally 
different from the assumption of "forced accretion of cold matter". 

Our 1-D stellar evolution models of the central remnant, i.e. 
the cold core and the hot envelope, which include the effects of 
rotation, indicate that there are three necessary conditions for the 



merger remnant to avoid off-center carbon ignition such that a SN 
la may be produced: 

(i) The local peak of temperature of the merger remnant at the 
interface between the core and the envelope must be lower than the 
critical temperature for carbon ignition (T p < Tc-i g ). 

(ii) The time scale for angular-momentum loss from the central 
remnant by must be larger than the neutrino cooling time scale at 
the interface (rj > t u ,p). 

(iii) Mass accretion from the Keplerian disc onto the cen- 
tral remnant must be sufficiently slow (M acc < 5 x 
10- 6 ...10~ 5 M yr" 1 ). 

Our new SPH simulation confirms that at least the first condition 
(T p < Tc-ig) should be fulfilled in the CO white dwarf binary 
considered. 

As emphasized in Sect. 13. 11 our 1-D models significantly un- 
derestimate the effect of the centrifugal force on the stellar structure 
in the rapidly rotating outermost layers. However, since thermal 
diffusion always dominates over both neutrino cooling and com- 
pressional heating in the outer envelope (p < 10 5 ...10 6 g cm -3 ) 
above the interface, as shown in Figs.[7]and[T2l the detailed struc- 
ture of the rapidly rotating outermost layers above the interface may 
not significantly affect our results on the thermal evolution of the 
merger remnant, as long as the angular momentum of the envelope 
is not lost faster than the local neutrino cooling time scale at the 
interface. On the other hand, mass accretion from the Keplerian 
disc should occur preferentially along the equatorial plane of the 
envelope. As shown in the SPH simulation, the envelope is more 
extended along the equatorial plane, where most angular momen- 
tum is deposited, than along the polar axis, and the resultant com- 
pressional heating must be much weakened, compared to the case 
of our 1-D models. The enhanced role of rotation must thus help 
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to increase the critical mass accretion rate for inducing off-center 
carbon ignition, in favor of producing a Type la supernova. 

We have concluded that the loss of angular momentum on a 
short time scale (rj < r u , P ~ 10 4 ...10 5 yr) may induce off-center 
carbon ignition even when T max ,init < Tc-i s . Rapidly rotating 
compact stars may experience loss of angular momentum by gravi- 
tational wave radiation, due to either the bar-mode instability or the 
r-mode instability. The onset of the dynamical or secular bar-mode 
instability requires a very high ratio of the rotational energy to the 
gravitational energy: E ro t/E grav > 0.2 for the dynamical bar- 
node instability, and E ro t/E &raY > 0.14 f or the secular bar-mode 
instability (e.g. Shapiro & Teukolskv 1983h . As both our 1-D mod- 
els and SPH simulation give a value of E rot / E grav that is much 
lower (about 0.06 - 0.07) than those critical limits, the bar-mode in- 
stability may not be relevant. The r-mode instabili ty may operate , 
in principle, even with suc h a low E YO t / E gra v (lAnderssonlll998l ; 
Friedman & Morsink 1998). However, we estimate that the growth 
time of the r-mod e instability (r r ), u sing our central remnant mod- 
els and following iLindbloml (L1999), is > 10 6 yr, which is much 
longer than the local neutrino cooling time scale (r u ,p ~ 10 4 yr). 
Alternatively, angular momentum might be transported from the 
accreting star into the Kepleria n disc when the accreto r reaches 
critical rotation. Calculations by ISaio & Nomotol (120041) indicate, 
however, that the decrease of the total angular momentum due to 
such an effect is not significant in accreting white dwarfs. In con- 
clusion, neither gravitational wave radiation nor outward angular- 
momentum transport is likely to lead to a rapid loss of angular 
momentum from the central remnant such that rj < r u ,p, unless 
magnetic torques are important. 

The central remnant may be enforced to rotate rigidly on a 
short time s cale in the presence of strong magnetic torques (cf. 
ISpruitll2002h . The central remnant in both our SPH simulation and 
1-D models has J to t > 10 50 erg s, which is significantly higher 
than the maximum limit a rigidly rotating white dwarf can retain, 
as shown in Fig.[l4] This means that if magnetic torques led to rigid 
rotation, a large amount of angular momentum should be trans- 
ported into the Keplerian disc (Case a in Fig. [14]), or mass shed- 
ding of super-critically spun-up layers should occur from the cen- 
tral remnant (Case b in Fig. [14]). In Case a, the local density around 
the interface should increase by several factors by the time when 
the central remnant reaches rigid rotation as implied by Fig. \\5\ 
Off-center carbon ignition might be inevitable in this case due to a 
resultant high effective accretion rate, if the time for angular mo- 
mentum redistribution were shorter than the local cooling time due 
to neutrino losses. In Case b, on the other hand, the local density 
at the interface might not increase if mass shedding from the cen- 
tral remnant occurred at a sufficiently high rate. Therefore, the role 
of magnetic fields in the merger evolution remains uncertain at the 
current stage and is a challenging subject for future work. 

The coalescence of more massive double CO white dwarf bi- 
naries is likely to result in a higher maximum temperature due to the 
enhanced role of gravity. Consequently, given the important role of 
the maximum temperature in the merger remnant for its final fate, 
less massive binary CO white dwarfs may be favored for the pro- 
duction of SNe la from such a channel. 

We note that there are number of potentially important factors 
that have not been included in either the present study or previous 
simulations. These include the following points: 

(i) The previous and present simulations assumed that white 

dwarfs are cold prior to the merging process. However, 

llben. Tutukov & Fedoroval dl998h point out that tidal interactions 




Figure 14. Evolution of the central remnant in the mass - angular momen- 
tum plane. The thick solid curve shows the angular momentum of a rigidly 
rotating white dwarf with critical rotation at the surface as a function of the 
white dwarf mass. The thick dashed curve and the thick dot-dashed curve 
give the critical angular momentum for a differentially rotating CO white 
dwarf to reach carbon ignition at the center (p c = 2 x 10 9 g cm -3 ), and 
electron-capt ure induced collapse (p c = 10xl0 10 g cm -3 ), respectively, 
according to lYoon &~L anger (2005). A SN la explosion is expected in the 
hatched region. The filled circle denotes the initial model of the central rem- 
nant in Seq. Sa. The evolution of the central remnant in Seqs Sa 8, Sa9 and 
Sal 1 is shown by the thin dotted curves, as indicated. The thin solid curves 
denote possible evolutionary paths of the central remnant with strong mag- 
netic torques that may enforce rigid rotation, with loss of angular momen- 
tum but without mass shedding (Case 'a'), and with both loss of angular 
momentum and mass shedding (Case 'b'). See the text for more details. 
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Figure 15. The density profile in the initial model of the central remnant 
in Seqs Sa (solid curve), and in a corresponding hot (T c = 10 8 K) white 
dwarf model that rotates rigidly at critical rotation at the surface (dashed 
curve). 



might heat up the white dwarfs as the orbit shrinks, which could 
weaken the gravitational potential of the primary. Furthermore, as 
the temperature of white dwarfs is a function of their age, younger 
progenitors should have more extended envelopes, which may re- 
sult in a lower T p . 

(ii) A thin hydrogen/helium envelope must be present initially in 
both the primary and the secondary. As hydrogen or helium should 
ignite at a much lower temperature than carbon, the influence of 
the release of nuclear energy during the merger process may be 
even more important than shown in the existing SPH simulations, 
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which is likely to lower T p . Furthermore, neutrino losses, which 
were neglected in the present study, would also tend to reduce T p . 

(iii) At a given total maSS (_/\/ftot — -^primary ~h -A^secondary )» 

different mass ratios of the white dwarf components (q = 
Mgecondary/Mprimary) must result in different merger structures. 

(iv) A lower q at a given M to t may not only lead to a stronger 
gravitational potential of the primary, but also to a lower mass 
accre tion rate during the dynamical mass transfer dGuerrero et al.1 
120041) . As the former and the latter will tend to increase and de- 
crease T p , respectively, quantitative studies are necessary to predict 
how T p will change with q. 

Finally, another important ingredient that needs to be consid- 
ered is thermal diffusion during the dynamical evolution. As shown 
above, the mass-accretion rate from the Keplerian disc onto the en- 
velope of the central remnant is one of the most important factors 
that critically determine the final fate of double CO white dwarf 
mergers. The accretion rates depend on the structure of the Kep- 
lerian disc at thermal equilibrium, which can be only understood 
by including thermal diffusion in future simulations. But here we 
emphasize again that the accretion rates from a centrifugally sup- 
ported Keplerian disc should be significantly lower than those from 
a pressure- supported thick disc that was previously assumed, which 
opens the possibility for at least some double CO white dwarf 
mergers to produce SNe la. 
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